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1 Introduction 



There are strong and well known indications that the three ordinary (active) neutrinos 
Ua (a = e, /i, r) are mixtures of different mass eigenstates uj, with j = 1, 2, 3, . . . (for 
a review see e.g. [|l|). Active neutrinos can mix among themselves, and possibly only 
among themselves, in which case j = 1, 2, 3. However, it is also possible that there are 
mixings of the active neutrinos with new sterile ones, Ug, that do not have any direct 
coupling to W and Z bosons mediating weak interactions. Sterile neutrinos have 
repeatedly been suggested as solutions of various anomalies in neutrino experiments. 
In particular, if all the present day data are correct, then at least one sterile neutrino 
must participate. For an explanation of the known neutrino anomalies, the mass 
differences should be very small, it should lie in the eV range or even in the sub-eV 
range. However, if we admit that there may be new neutrino states mixed to the 
usual ones, we should keep an open mind to possible values of their masses. It is an 
interesting question which values of masses and mixing angles can be excluded by 
direct experiments as well as by cosmology and astrophysics. The last two topics are 
addressed in the present paper which was partly stimulated by the recent limits on 
mixing of Ur with a heavy sterile neutrino found by the NOMAD collaboration 0, 
and partly by our paper P] where we have found that cosmology and astrophysics 
practically exclude the interpretation of the KARMEN anomaly by a 33.9 MeV 
neutrino mixed with z/^ (astrophysical and cosmological limits on 33.9 MeV sterile 
neutrino were also considered in ref. ||5|). According to a statement of the KARMEN 
collaboration made at Neutrino 2000 the anomaly was not observed in upgraded 
detector KARMEN 2, but the question still remains which area in the mass-mixing- 
plane can be excluded. 

In what follows we will derive the bounds on masses and mixings that follow from 
big bang nucleosynthesis (BBN) and from the duration of the supernova (SN) 1987A 
neutrino burst. We essentially repeat our previous work 0, but now lifting the restric- 
tion on a specific value of mass and lifetime of the heavy neutrinos and considering 
mixings with all active states Ug, and Vr- We assume that the heavy neutrinos 
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have Dirac mass. In the case of Majorana mass the hmits would be shghtly weaker. 



2 Preliminaries 

We assume that the sterile neutrino mixes predominantly with only one active flavour 
'^a = ^'/i or ^T- The mixed flavour states are expressed in terms of the mass 
eigenstates and the mixing angle 6 as 

Pa = cos 9 ui + sin 9 U2 , 

Vs = — sin 9 ui + cos 9 1^2 , (1) 

where ui and z/2 are assumed to be the light and heavy mass eigenstates, respectively. 
This mixing couples the heavier neutrino to the Z", allowing for the decay channel 

iy2^^i+i + i, (2) 

where ui is mostly an active flavour and i is any lepton with a mass smaller than the 
mass m2 of the heavy neutrino. We assume that 1712 < 2m^ so that the decay into 
flfi and TT is kinematically forbidden. If the active neutrino mixed with Us is either 
Ufj^ or then we can translate the mixing angle into lifetime as 

(Ms/lO MeV)5 sm^29 ' ^ ' 

For the mixing with the numerator is 0.7 sec; it is is different because of the 
presence of charged-current interactions. 

A sterile neutrino mixed with z/^ could be observed in direct experiments, in par- 
ticular where upper bounds on neutrino masses are obtained (see the list of references 
in [0). The most accurate limits of course exist for ^ roughly m^^ < 3 eV. 
However, these experiments cannot help much in eliminating a heavy sterile neutrino 
because they are not sensitive to the mass range > 10 MeV which we consider. 
Such heavy neutrinos are simply not produced in beta-decays and their impact is only 
indirect, e.g. they renormalize vector and axial coupling constants. 

The upper limit on the tau-neutrino mass obtained by the ALEPH Collabora- 
tion could be translated into limits on mass/mixing of Ur with Ug. However, one 



would need to reanalyse the data under the assumption of two (or several) mixed 
neutrinos, taking account of the sensitivity to measure the energy spectrum for dif- 



ferent values of m^. The bound obtained in Ref. [jTT| based on the assumption that 
the average neutrino mass (m) = cos^ 6 mi + sin^ 9 m2 should be smaller than the ex- 
perimental upper limit is too naive and can serve only as a rough order-of-magnitude 
estimate. Moreover the NOMAD bounds are much more restrictive and thus we 
will not pursue the subject of the ALEPH bounds here. 

3 Cosmological production and freeze-out of heavy 
sterile neutrinos 

In the early universe sterile neutrinos are produced through their mixing with the 
active ones. The production rate for relativistic Vg (i.e. for > 1712) can be approxi- 
mately estimated as [|l^]: 



where H is the Hubble expansion parameter, Ty is the plasma temperature equal 
to the photon temperature. is the decoupling temperature of active neutrinos, 
taken to be 3 MeV, and 6'm is the the mixing angle in the medium. According to the 



calculations of Ref. [|13[ one finds in the limit of small mixing: 



. „^ sin 26* 

sin 20i\/r ~ f5) 

1 + 0.76 X 10-19 (r^/MeV)6(5mVMeV')-i ' ^' 

The matter effects become essential for 

Ty > 1.5 X 10^MeV(W/MeV2)^/^ (6) 

For Tg/H > 1, sterile neutrinos were abundantly produced and their number 
density was equal to that of light active neutrinos, at least during some epoch. The 
production rate reaches a maximum at Tniax = 1.28 X m\6m^ /MeV^y/^ MeV. For 
the masses that are considered below, T^ax is well above the neutrino mass. 

If the equilibrium number density of sterile neutrinos is reached, it would be 
maintained until Ty ^ 4(sin2^)-^/^ MeV. This result does not depend on the heavy 



neutrino mass because they annihilate with massless active ones, 1^2 + J^a ^ M- The 
heavy neutrinos would be relativistic at decoupling and their number density would 
not be Boltzmann suppressed if, say, Tf > Ms/2. This gives 

sm^2e{6m^/MeVy^^ < 500 . (7) 

If this condition is not fulfilled the impact of Us on BBN would be strongly diminished. 
On the other hand, for a sufficiently large mass and non-negligible mixing, the z/2 
lifetime given in Eq. (^) would be quite short, so that they would all decay prior to 
the BBN epoch. (To be more exact, their number density would not be frozen, but 
follow the equilibrium form oc e~*^^/^^.) 

Another possible effect that could diminish the impact of heavy neutrinos on BBN 
is entropy dilution. If z/2 were decoupled while being relativistic, their number density 
would not be suppressed relative to light active neutrinos. However, if the decoupling 
temperature is higher than, say, 50 MeV pions and muons were still abundant in the 
cosmic plasma and their subsequent annihilation would diminish the relative number 
density of heavy neutrinos. If the decoupling temperature is below the QCD phase 
transition the dilution factor is at most 17.25/10.75 = 1.6. Above the QCD phase 
transition the number of degrees of freedom in the cosmic plasma is much larger 
and the dilution factor is approximately 5.5. However, these effects are essential for 
very weak mixing, for example the decoupling temperature would exceed 200 MeV if 
sin^ 2^^ < 8 X 10~^. For such a small mixing the life-time of the heavy z/2 would exceed 
the nucleosynthesis time and they would be dangerous for BBN even if their number 
density is 5 times diluted. 

Sterile neutrinos would never be abundant in the universe if Tg/H < 1. In fact 
we can impose a stronger condition demanding that the energy density of heavy 
neutrinos should be smaller than the energy density of one light neutrino species at 
BBN (T ~ 1 MeV). Taking into account a possible entropy dilution by factor 5 we 
obtain the bound: 

(SmyUeV^) sin^ 26 < 2.3 x 10"^ . (8) 
Parameters satisfying this conditions cannot be excluded by BBN. A more detailed 
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consideration permits one to impose a somewhat better limit, but we will not go into 
such detail here. 

Keeping these restrictions in mind we proceed to derive bounds on Mg and the 
mixing angle 9 demanding that the influence of heavy neutrinos on BBN should not 
be too strong. 

4 Big-Bang Nucleosynthesis 

In general, if a neutrino has a mass exceeding an MeV then it may influence the light 
element abundances unless it decays much before the onset of nucleosynthesis (see 



e.g. [14, 15 1 and references therein). A heavy unstable sterile neutrino will affect BBN 
in several ways. The main effect is through the increased energy density which leads to 
a faster expansion and hence an earlier freeze-out of the n/p-ratio. However, also the 
decay products must be taken into account since the electron neutrinos directly enter 
the n-p reactions. Moreover, if the decays into the e+e^ channel the temperature 
evolution is altered (see 0, |l^ for discussion). 

The calculations are described in detail in Ref. 0; here we only briefly describe 
the basic steps. First we introduce the dimensionless variables 

X = rriQait) and y=pa(t) , (9) 

where a{t) is the cosmic scale factor and mo is an arbitrary normalisation factor that 
we have chosen as mo = 1 MeV. We also introduce a dimensionless temperature 
T = T^/tuq and measure all masses in units of mo. In terms of these variables the 
Boltzmann equations describing the evolution of the sterile neutrino, Ug, and the 
active neutrinos, z/^, can be written as 



fus - fusj 1.48 X / 10.75' 



1/2 



r.,/sec {Txy \g.{T) J 

Mg ^ 3 X 2^r3 / 3C(3) ^ Ttt^ / ^ plj 



(10) 



M3 V 4 144 V M2 ?>E,M^J^ 
dJU^ = Da{x,y,Ms) + Sa{x,y), (11) 
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where f^^ = (e^/^ + iy\ 5^ = U - (e^ + 1)^', and E,^ = ^M^ + {y/xf and 
Vus = y/x are the energy and momentum of z/^. Further, (?=„(T) = ptot/(7r^T^/30) is 
the effective number of massless species in the plasma determined as the ratio of the 
total energy density to the equihbrium energy density of one bosonic species with 
temperature T^. 

The scattering term in Eq. (0) comes from interactions of active neutrinos be- 
tween themselves and from their interaction with electrons and positrons. For Ue it 
has the form 

/in 75\^/^ 

S^^{x,y) =0.26(^^j {l + 9l + gl){y/x') 

X J dy2yl6fuXx,y2) + ^ / dy2yl [5fu^{x,y2) + 5/^, (a;, 1/2)) 

where = sin^ 6w + \, while qr = sin^ 9\y- The corresponding term for z/^ and is 
given by Eq. (|1^) with the exchange Ql Ql = sin^ ^'ly — |. 

The decay term in Eq. (|ll]) comes from the decay of heavy sterile neutrino ac- 
cording to reactions Eq. (^. This term depends on the mixing channel. For Uj. ^ 
mixings we have 




X ("^. - "2) iM.xl2 - y) , (14) 

where n^^{x) is the number density of Ug and 6{y) is the step function which ensures 
energy conservation in the decay. For z/^ ^ Ug mixing these terms come from Eq. (|13D 
and Eq. (|^) by exchange of ^ Ur- 
For Ue <-> Us mixing those terms are 
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X {n,^-nl^)9iMj2-y/x) , (16) 

More details on the derivation of Eqs. (|TU|,|TT]) can be found in Ref. 

The temperature evolution is governed by the equation for energy conservation: 

x^ = -3ip + p), (17) 

where p and p are energy density and pressure. 

For the solution of Eqs. (|10],|lT],|l7|) we divide the time region into 3 parts. At early 
times, when the light neutrinos are kept in tight equilibrium with the electromagnetic 
plasma, we need only to solve the equation for the heavy sterile neutrinos, Eq. ([T0|). 
At intermediate temperatures we solve the full set of equations, taking into account 
the changed temperature evolution and the non-equilibrium distribution functions 
of the neutrinos. At late times, long after the neutrinos have decoupled and the 
heavy neutrinos have disappeared, we solve only the kinetic equations governing the 
n-p reactions needed for the nucleosynthesis code. We use a 800 point grid in the 
momentum region big enough to encompass the decay products {0 < y < M<,x/2), 
and for the BBN calculations we use ?7io = 5. 

The results of the calculations have been imported into the modified Kawano nu- 
cleosynthesis code [|T7|, and the abundances of all light elements have been calculated. 
At each time step x we find the corresponding photon temperature and total energy 
density. Furthermore we integrate the kinetic equation governing the n/p evolution 
taking into account the distorted spectrum of z/g. 

For a given mass, Mg, we calculate the helium abundance, Y, for various mixing 
angles (or lifetimes) and use the approximate formula AY = 0.013 AA^ to find the 
minimum allowed mixing angle (Fig. |l|) or the maximum allowed lifetime (Fig. 0) for 
that mass. We do this for both AA^ = 0.2 and 1.0. On the same figures we plot 
guide-the-eye lines, that are fitted in the form 

(sin2^)^i, = siM; + S2 (18) 
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Figure 1: Minimum mixing angle between sterile and active neutrinos, allowed by 
BBN, as a function of the heavy neutrino mass, both for an optimistic bound, AN = 
0.2, and for a conservative bound, AN = 1.0. The upper panel corresponds to Ug 
or <-> Vs mixings, while the lower one to ^ J^s mixing. 
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Figure 2: Maximum lifetime of sterile neutrino, allowed by BBN, as a function of the 
heavy neutrino mass, both for an optimistic bound, AA^ = 0.2, and for a conservative 
bound, AN = 1.0. The upper panel corresponds to Us or Us mixings, 

while the lower one to ^ J^s mixing. 
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for the minimum allowed mixing and 

W = h Mf + t2 (19) 



for the maximum allowed life-time Tmax- Fitting parameters in Eqs. (|TB|,|IUp are pre- 
sented at the Table |l]. 
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(sin^6')min fitting parameters. 


Tmax fitting parameters. 


Sl 




a 
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-3.727 
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0.2 
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1218 


0.0513 


-2.658 


Vs - 


1 


66.33 


-1.05 - 10-5 


-3.070 


1699 


0.0544 


-2.652 



Table 1: Fitting parameters for minimum allowed mixing (sin^ ^)min fitted by Eq.(|l^) 
and maximum allowed life-time Tmax fitted by Eq.([l9|). 



For small masses we followed Ref. and used the expansion in a small parameter, 
A = T - X — 1, to describe the temperature evolution, whereas for bigger masses 
we used a more correct, and somewhat more complicated, treatment of the photon 
temperature from Eq. (PTf). The results are in perfect agreement for Ms = 60 MeV, 
and for Mg = 140 MeV we find the maximal allowed lifetime about 0.05 sec, in fair 
agreement with the approximate fitting formulae. 

It is important to keep in mind, that the bounds obtained above are rather conser- 
vative, since all the approximations used in the derivation lead to slightly weakened 
bounds, as described in Ref. [0]. 
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5 Decay 1^2 + tt^ 

The results obtained in the previous section are vahd when the mass of the heavy 
neutrino is smaller than 140 MeV. For such low masses the dominant decay modes 
of i>2 would be into electrons and light neutrinos. A possible decay mode including 
muons which is possible for Ms > 105 MeV is suppressed if Mg < 140 MeV. However, 
for Ms > 135 MeV the decay channel 

z/2 ^ 7r° + Ua (20) 

becomes open and strongly dominating. The life-time of vr'' vv was calculated 
in refs. [p!8| , p!9|] . We can translate their results for the decay ( pOD and find for the 



life-time 



r 



GlMs{Ml-ml)fl^xYv^e 



5.8- 10"^ sec 



sin^^ 



MsiMl - m 



2 ™2n"I-1 



ml 



(21) 



167r 

where = 135 MeV is the vr^-mass and fn = 131 MeV is the coupling constant for 
the decay 7r+ ^ yU + z/^. 

Immediately after Mg becomes bigger than the two-body decay becomes the 
main one and all other channels can be neglected. We cannot directly apply our 
numerical program (that was made for a 33.9-MeV neutrino) to this case, so we will 
instead make some simple order of magnitude estimates for the impact of very heavy 
z/2 on BBN. One can roughly conclude that for the life-time of z/2 smaller than 0.1 sec, 
and corresponding cosmological temperature higher than 3 MeV, the decay products 
would quickly thermalize and their impact on BBN would be small. For a life-time 
larger than 0.1 sec, and T < 3 MeV, we assume that thermalization of neutrinos 
is negligible and approximately evaluate their impact on BBN. If Ug is mixed with 
z/^ or Uj. then electronic neutrinos are not produced in the decay (pOD and only the 
contribution of the decay products into the total energy density is essential. As we 
have already mentioned, non-equilibrium z/g produced by the decay would directly 
change the frozen ra/p-ratio. This case is more complicated and demands a more 
refined treatment than what is presented in this section. 
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The 7r° produced in the decay (|20|) immediately decays into two photons and they 
heat up the electromagnetic part of the plasma, while neutrinos by assumption are de- 
coupled from it. We estimate the fraction of energy delivered into the electromagnetic 
and neutrino components of the cosmic plasma in the instant decay approximation. 
Let Ts = Us/rio be the ratio of the number densities of the heavy neutrinos with 
respect to the equilibrium light ones, uq = 0.09T^. The total energy of photons and 
e"'"e~-pairs including the photons produced by the decay is 

P^m = --T + (22) 

while the energy density of neutrinos is 

The effective number of neutrino species at BBN can be defined as 

= (24) 

Because of the stronger heating of the electromagnetic component of the plasma by 
the decay products, the relative role of neutrinos diminishes and Kjf-^^^ becomes 
considerably smaller than 3. If are decoupled while relativistic their fractional 
number would be = 4 (two spin states and antiparticles are included). Possible 
entropy dilution could diminish it to slightly below 1. Assuming that the decoupling 
temperature is = 3 MeV we find that Kl'^ff^ = 0.6 for Ms = 150 MeV and 
K^ff^ = 1.3 for Ms = 200 MeV if the frozen number density of Us is not diluted 
by the later entropy release and remains equal to 4. If it was diluted down to 1, 
then the numbers would respectively change to K^^^^^ = 1.15 for Mg = 150 MeV and 
ir^e//) ^ l Y for Ms = 200 MeV, instead of the standard K^^^^'^ = 3. Thus a very 
heavy Us would result in under-production of ^ife. There could, however, be some 
other effects acting in the opposite direction. 

Since Uf. decouples from electrons/positrons at smaller temperature than z/^ and 
Ut-, the former may have enough time to thermalize. In this case the temperatures of 
z/g and photons would be the same (before e"'"e~-annihilation) and the results obtained 
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above would be directly applicable. However, if thermalization between z/g and was 
not efficient, then the temperature of electronic neutrinos at BBN would be smaller 
than in the standard model. The deficit of Vf, would produce an opposite effect, 
namely enlarging the production of primordial ^ife, because it results in an increase 
of the ra/p-freezing temperature. This effect significantly dominates the decrease of 
K^^^"^ discussed above. Moreover even in the case of the decay z/2 ^ 7r° + v^^r-, when 
Vf, are not directly created through the decay, the spectrum of the latter may be 
distorted at the high energy tail by the interactions with non-equilibrium and Vy, 
produced by the decay. This would result in a further increase of ^iJe-pro duct ion. In 
the case of direct production of non-equilibrium z/g through the decay z/2 ^ tt*^ + z/g 
their impact on n/p ratio would be even much stronger. 

To summarise, there are several different effects from the decay of Vg into 7r° 
and V on BBN. Depending upon the decay life-time and the channel these effects 
may operate in opposite directions. If the life-time of V2 is larger than 0.1 sec but 
smaller than 0.2 sec, so that and z/g establish equilibrium the production of '^He is 
considerably diminished so that this life-time interval would be mostly forbidden. For 
life-times larger than 0.2 sec the dominant effect is a decrease of the energy density 
of z/g and this results in a strong increase of the mass fraction of ^ife. Thus large life- 
times should also be forbidden. Of course there is a small part of the parameter space 
where both effects cancel each other and this interval of mass/mixing is allowed. It 
is, however, difficult to establish its precise position with the approximate arguments 
described above. It would be a matter of separate and rather complicated non- 
equilibrium calculations. 

Thus, in the case of z/^ ^ z/^^^ mixing for Mg > 140 MeV we can exclude the 
life-times of Ug roughly larger than 0.1 sec, except for a small region near 0.2 sec 
where two opposite effects cancel and the BBN results remain undisturbed despite 
the presence of sterile neutrinos. Translating these results into mixing angle according 
to Eq. (piD , we conclude that mixing angles sin^ 9 < 5.8 ■ 10^^m.,^/Ms/ [[Ma/m.,^Y — 1) 
are excluded by BBN. Combining this result with Eq.(|^) we obtain the exclusion 
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region for > lAOMeV: 

o MeV^ 9 „ o 1 , , 

51 aO"^ — -— < sin^^ < 5.8- 10"* -. 25 

In the case of Vs ^ z/g mixing for Ms > lAOMeV the hmits are possibly stronger, 
but it is more difficult to obtain reliable estimates because of a strong influence of 
non-equilibrium z/g produced by the decay on neutron-proton reactions, and we will 
postpone this problem for future investigations. 

6 Supernova 1987A 

In Ref. we have explained that massive sterile neutrinos can be constrained by the 
duration of the SN 1987A neutrino burst. If the sterile states live long enough and 
interact weakly enough to escape from the SN core, the usual energy-loss argument 
constrains the allowed interaction strength (in our case the active-sterile mixing an- 
gle). If the sterile states are so short lived or so strongly interacting that they decay 
or scatter before leaving the SN core, they contribute to the transfer of energy and 
thus still accelerate the cooling speed of the SN core. 

The arguments given in Ref. |^ pertain directly to the present long as the 

sterile neutrino mass is not too large for its production to be suppressed in the relevant 
region of a SN core. Taking a typical temperature to exceed 30 MeV, the average 
thermal neutrino energy exceeds about 100 MeV so that threshold effects should 
not become important until the mass significantly exceeds this limit. Therefore, we 
believe that the SN limits pertain at least to masses up 100 MeV. For such a mass and 
maximal mixing, the lifetime is about 10 /iS, for relativistic particles corresponding 
to a distance of about 3 km, i.e. of order the SN core radius, and much larger than 
the standard mean free path of active neutrinos of a few meters. Therefore, for the 
entire range of interesting masses and mixing angles the lifetime far exceeds what is 
necessary for the sterile neutrinos to transfer energy within the SN core in the spirit 
of the energy-transfer argument. 

In summary, the arguments of Ref. imply that for sterile neutrino masses below 
about 100 MeV the approximate range 3 x 10~* < sin^(2^) < 0.1 is excluded. 
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7 Conclusion and Summary 



In this work we have found cosmological and astrophysical bounds on possible mixings 
of a heavy sterile neutrino with mass 10 MeV < Ms < 140 MeV with any active 
flavour z/e, or z/^. Our results are summarised in Fig. |^. The region between the 
two horizontal lines running up to 100 MeV are excluded by the duration of the 
neutrino burst from SN 1987A. A more detailed discussion would probably permit 
one to expand this region both in the horizontal and vertical directions. 

The BBN limits are presented by the two upper dashed curves, both corresponding 
to the conservative bound of one permitted extra neutrino species. The curve for 
v^^r mixing is slightly higher than the curve for mixing. The lower dashed curve 
describes our approximate estimate of the efficiency of sterile neutrino production in 
the early universe Eq. @. Below this curve the heavy z/2 are very weakly produced 
and their impact on BBN is not essential. Let us note that relation Eq. is fulfilled 
for the entire BBN-excluded region of Fig. ^. In other words, sterile neutrinos are 
relativistic at decoupling and the Boltzmann suppression factor is not essential. 

We have also made an order of magnitude estimates of the BBN-influence of 
heavier neutrinos with masses 140 MeV < Ms < 200 MeV for Vg v^^t mixing 
channel. We have found, that mixing angles in the range given by Eq. (|25|) are excluded 
by BBN. We do not present these order of magnitude estimates in our Fig. ^ since 
they are less precise than the limits presented for Mg < 140 MeV. We believe that the 
BBN bounds also could be improved considerably with more detailed calculations. 
For the ^ mixing a part of the sin^ ^-Mg-plane is excluded by the NOMAD 
experiment 0. Their bound is also presented in Fig. ^ 

We have done the calculations for the Dirac neutrinos. For the Majorana case 
the cosmological number density of the heavy neutrinos would be twice smaller and 
corresponding bounds would be changed in an evident way and rather weakly. 

It is noteworthy that the impact of sterile neutrinos on BBN is very sensitive to 
the deviation from thermal equilibrium of both sterile neutrinos as well as their de- 
cay products, especially into the i/e-channel. In the usual approximate calculations. 
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Figure 3: Summary of our exclusion regions in the sin^ 6'-Ms-plane. SN 1987A ex- 
cludes all mixing angles between two solid horizontal lines. BBN excludes the area 
below the two upper dashed lines if the heavy neutrinos were abundant in the early 
universe. These two upper dashed lines both correspond to the conservative limit 
of one extra light neutrino species permitted by the primordial ^He-abundance. The 
higher of the two is for z/^^^ mixing, and the slightly lower curve is for z/g mixing. 
In the region below the lowest dashed curve the heavy neutrinos are not efficiently 
produced in the early universe and their impact on BBN is weak. For comparison we 
have also presented here the region excluded by NOMAD Collaboration for the 
case of i^s ^ Vt mixing. 



kinetic equilibrium is assumed. This assumption very much simplifies the calcula- 
tions because the integro-differential kinetic equations are reduced to simple ordinary 
differential equations. However, we have found that this approximation deviates 
substantially from the more exact treatment presented here. Even more accurate cal- 



culations are feasible with the technique developed in our previous works [T^, [16 



where an accuracy at the sub-percent level was achieved, but this is a difficult and 
time-consuming problem. We would like to postpone it until the existence of heavy 
sterile neutrinos becomes more plausible. 
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